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Ground state selection in XY pyrochlore antiferromagnets R 2 Ge 2 0 7 (R = Er, Yb) 


Z. L. Dun, 1 X. Li, 2 R. S. Freitas, 3 E. Arrighi, 3 C. R. Dela Cruz, 4 M. Lee, 5 ’ 6 E. S. Choi, 6 
H. B. Cao, 4 H. J. Silverstein, 7 C. R. Wiebe, 6 ’ 7 ’ 8 ’ 9 J. G. Cheng, 2 and H. D. Zhou 1 - 6 

department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996-1200, USA 
2 Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
3 Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 05314-970, Sao Paulo, SP, Brazil 
* Quantum Condensed Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37381, USA 
5 Department of Physics, Florida State University, Tallahassee, FL 32306-3016, USA 
6 National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL 32310-3706, USA 
7 Department of Chemistry, University of Manitoba, Winnipeg, MB R3T 2N2 Canada 
8 Department of Chemistry, University of Winnipeg, Winnipeg, MB, R3B 2E9 Canada 
9 Canadian Institute for Advanced Research, Toronto, Ontario, M5G 1Z7, Canada 

(Dated: August 21, 2015) 

Elastic neutron scattering, ac susceptibility, and specific heat experiments on the pyrochlores 
Er 2 Ge 207 and Yb 2 Ge 207 show that both systems are antiferromagnetically ordered in the r 5 
manifold. The ground state is a ^3 phase for the Er sample and a ^2 or ^3 phase for the Yb sample, 
which suggests “Order by Disorder” (ObD) physics. Furthermore, we unify the various magnetic 
ground states of all known R 2 B 2 O 7 (R = Er, Yb, B = Sn, Ti, Ge) compounds through the enlarged 
XY type exchange interaction J± under chemical pressure. The mechanism for this evolution is 
discussed in terms of the phase diagram proposed in the theoretical study [Wong et al., Phys. Rev. 

B 88, 144402, (2013)]. 

PACS numbers: 75.10.Jm, 61.05.fm, 75.40.-s 


The pyrochlores R2B2O7 (R: rare earth elements, B: 
transition metals) have been a hot topic due to their 
emergent physical properties based on the geometrically 
frustrated lattice mm- Recent interest in pyrochlores 
is focused on systems with effective spin-1/2 R 3+ ions 
012 , in which the crystal electric field (CEF) normally 
introduces a well-isolated Kramers doublet ground state 
with easy XY planar anisotropy Ei- In these XY 
pyrochlores, the anisotropic nearest neighbor exchange 
interaction J ex = (J zz , «/±, Jz±, J±±) between the R 3+ 
ions, plus the strong quantum spin fluctuations of the 
effective spin-1/2 moment, stabilize various exotic mag¬ 
netic ground states [3j- 

Er2Ti207 and Yb2Ti2©7 are two celebrated examples 
of the effective spin-1/2 XY pyrochlores. For Yb2Ti207, 
the local [111] Ising-like exchange interaction J zz is con¬ 
siderably larger than the XY planar interaction J± [7j. 
An unconventional first order transition is observed [8], 
which has been proposed to be a splayed-ferromagnet 
(SF) state with Yb 3+ spins pointing along one of the 
global major axes with a canting angle [9j. Slight dis¬ 
order between the Yb and Ti sites leads to a possible 
quantum spin liquid state [TUI HI] . For Er 2 Ti 2 07, the 
Er 3 + spins are energetically favored to lie within the 
local XY plane due to the dominating J±, in which a 
continuous U(l) symmetry is preserved in the Hamilto¬ 
nian that allows the Er 3+ spins to rotate continuously 
in the XY plane [HHH]- Recently, both experimental 
and theoretical studies suggest that the quantum spin 
fluctuations lift the U(l) degeneracy with a small gap 
opening in the spin-wave spectrum and select an anti¬ 


ferromagnetic (AFM) ordering state (7/2) as the ground 
state for E^T^C^. This is the so called “order by dis¬ 
order” (ObD) mechanism [T414T8] . in which the ground 
state is selected through entropic effects. Meanwhile, an 
alternative CEF-induced energetic selection mechanism 
is proposed that will likewise result in the 7/2 state with 
similar value of the gap|19, |20j. 

These delicate magnetic ground states are fragile and 
easily affected by perturbations, such as chemical pres¬ 
sure. By replacing the Ti 4+ sites with the nonmag¬ 
netic Sn 4+ and Ge 4+ ions, the lattice parameter varies 
to changes the exchange interactions. As listed in Ta¬ 
ble I, for both Er2B2 07 and Yb2B207 series, the Curie 
temperature and ordering temperature increase with de¬ 
creasing lattice parameter. Moreover, their magnetic 
ground states are markedly different. Er2Sn207 does not 
show any long-range magnetic ordering down to 50 mK 
[21] but displays a spin freezing below 200 mK with the 
AFM Palmer-Chalker (PC) correlations [22]. It’s pro¬ 
posed that Er2Sn207 is approaching the 7/2/PC phase 
boundary where the selection of either state is weak [221 - 
[24]. Er2Ge2C>7 shows an AFM ordering[25] that is similar 
to Er 2 Ti 2 07. While a similar SF phase is observed for 
both Yb2Ti 2 07 and Yb2Sn 2 07 j2Glf28l . Yb2Ge207 strik¬ 
ingly displays AFM ordering at TV = 0.61 K [29] . So 
far, the exact nature of the magnetic ground states of 
Er2Ge207 and Yb2Ge2C>7 are not clear. Are they also 
selected by ObD mechanism? More importantly, while 
the theoretical studies [SHU HU have made significant 
efforts to unify the magnetic properties of Yb and Er- 
XY pyrochlores, unified magnetic phase diagrams have 
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TABLE I: Comparison between Er 2 B 2 07 and Yb 2 B 2 07 - 




Er 2 B 2 07 



Yb 2 B 2 07 

B site ion 

Sn 

Ti 

Ge 

Sn 

Ti 

Ge 

ir(b 4+ )(A) 

0.69 

0.605 

0.53 

0.69 

0.605 

0.53 

a( A) 

10.35 

10.07 

9.88 

10.28 

10.03 

9.83 

e C w{ K) 

-14 

-15.9 

-21.9 

0.53 

0.75 

0.9 

t n 


1.17 

1.41 

0.15 

0.24 

0.62 

Order type 

~(AFM) AFM 

AFM 

FM 

FM 

AFM 

Reference 

m 


m 

m 

m 

m 

Spin state 

-(PC) 


^3 

SF 

SF 

l±(or3) 

Reference 

m 

[15| this work 

m 

m 

this work 


not been experimentally achieved. 

In this letter, we studied the polycrystalline py- 
rochlores Er 2 Ge 207 and Yb 2 Ge 207 using elastic neutron 
scattering under magnetic fields, ac susceptibility, and 
specific heat measurements. We identified a ^3 phase for 
the Er sample and a ^2 or ^3 phase for the Yb sample 
(see Fig. 1 (e)(f) for their spin configurations), which 
suggest ObD mechanism. Furthermore, we unified the 
various magnetic ground states of all studied R 2 B 2 O 7 (R 
= Er, Yb, B = Sn, Ti, Ge) through the enlarged XY 
type exchange interaction J± under chemical pressure. 
We discussed this general rule in terms of the phase dia¬ 
gram proposed by Wong et al [24] . 

Experimental details are listed in the supplemental ma¬ 
terials. By comparing the neutron diffraction patterns 
measured at 3 K and 0.3 K (Fig. l(a, b)) for Er 2 Ge 207 , 
several magnetic Bragg peaks, such as ( 111 ) (220) (311), 
etc., are clearly observed at 0.3 K ( < T/v = 1.41 K). 
The refinements using the XY type AFM spin structure 
in the manifold, either ^2 or vp 3 (Fig. l(e, f)), fit these 
magnetic Bragg peaks well with a magnetic moment of 
3.23(6) fiB • In fact, all magnetic phases within the T 5 
manifold result in the same diffraction pattern and it’s 
impossible to distinguish them in powder samples with 
zero-field data. Fig .1(d) shows the field dependence of 
the (220) and (311) Bragg peaks’ intensities. The details 
are: (i) with H < 0.15 T, a magnetic domain alignment 
results in a quick drop of the ( 220 ) peak intensity with 
increasing field; (ii) between 0.15 and 2 T, the spins grad¬ 
ually rotates with the magnetic field but keeps the AFM 
nature; (iii) around a critical field H c = 2 T, the (220) 
Bragg peak’s intensity abruptly drops to a background 
value while the (311) Bragg peak’s intensity continuously 
increases. This demonstrates that above H c , Er 2 Ge 207 
enters a spin polarized state. The observed FM (400) and 
AFM (200) Bragg peaks on the pattern measured at H = 
5 T (Fig. 1(c)) suggest that this polarized state is similar 
to the SF state in the T 9 manifold. The refinement by 
assuming one single SF structure with the magnetic field 
applied along the global z axis (Fig .1(g))actually fits the 
powder average 5 T data well with the Er 3+ moment as 
M=(±1.42(2), ±1.42(2), 4.40(1)) /ip in the global coor¬ 
dinate frame. The double peak feature of the reported ac 
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FIG. 1: (color online) Elastic neutron scattering patterns and 
Rietveld refinements for Er 2 Ge 2 07 at (a)T = 3 K and H = 0 
T, (b)T — 0.3 K and II = 0 T, and (c)T — 0.3 K and H = 
5 T. (d) The field dependence of the (200) and (311) Bragg 
Peaks’ intensities measured at T = 0.3 K, the critical field H c 
is marked as the dash line. The spin configurations for (e)^ 2 , 
(f) , 03 and (g)splayed-ferromagnetic (SF) phases in the local 
coordination. 


susceptiblity data for Er 2 Ge 207 also confirmed the mag¬ 
netic domain alignment around 0.15 T and the critical 
field around 2 T [25] . 

It has been pointed out m that (i)for both ^2 and 
ips states, a multi-domain state with equal fraction of 6 
magnetic domains (plotted in the supporting material) at 
zero field will be expected, which give different intensities 
of the (220) Bragg peak; (ii)with the applied magnetic 
field in [ 110 ] direction, two domains with larger intensity 
will be selected if the ^2 phase is present HU EH- This 
will result in a ( 220 ) peak’s intensity jump, which has 
been exactly observed for Er 2 Ti 207 in the single crys¬ 
tal neutron diffraction experiments mm-, (iii) similarly, 
if the ip 3 state is selected, a decrease is expected for 
the ( 220 ) peak’s intensity since the two domains with 
lower intensities will be selected. In our neutron pow¬ 
der diffraction experiment by using a pelleted sample, 
the magnetic field was applied vertically such that it is 
perpendicular to the scattering plane. Then a similar se¬ 
lection rule would be expected in addition to a powder 
averaging effect (see detailed analysis in the supporting 
material). As shown in Fig. 1(d), the (220) peak’s in¬ 
tensity drops dramatically from 400 at 0 T to 250 counts 
at 0.15 T. This result suggests that Er 2 Ge 207 orders in 
the -03 phase. However, in order to provide unambiguous 
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FIG. 2: (color online) (a)Elastic neutron scattering pattern 
and Rietvelt refinement for Yb 2 Ge 207 at T = 0.3 K and H 
= 0 T. (b) The difference between the patterns measured at 
0.3 K (with H = 0 and 2 T) and 1.6 K. (c) The field depen¬ 
dence of the (200), (311) and (400) Bragg Peaks’ intensities 
at 0.3 K. (d) The ac susceptibility of Yb 2 Ge 207 at different 
temperatures. Insert: the dc magnetization measured at 0.6 
K. (e) The magnetic phase diagram of Yb 2 Ge 2 07 . 

evidences for the ^3 state, polarized neutron experiments 
on a single crystal sample are needed. 

Fig. 2(a) shows the neutron diffraction pattern mea¬ 
sured at 0.3 K (< T n = 0.62 K) for Yb 2 Ge 2 07 . Due to 
the small magnetic moment of the Yb 3+ ions, the mag¬ 
netic Bragg peaks are weak (as shown in the insert). The 
difference between the 0.3 K and 1.6 K patterns (Fig. 
2(b)) more clearly shows that the observed magnetic 
Bragg peaks’ positions and intensity ratios are very sim¬ 
ilar to those of Er 2 Ge 2 07 , which identifies Yb 2 Ge 2 07 5 s 
ground state as either ^2 or ^3 in the manifold. Re¬ 
finements based on these two spin structures give the 
same Yb 3+ moment of 1.06(7) /i#, which is consistent 
with the previous report (Yb 3+ ~ 1.15 /ig) [32 . 

With an applied magnetic field on Yb 2 Ge 2 07 (Fig. 
2 (c), the ( 220 ) peak’s intensity decreases quickly around 
0.2 T, which indicates a critical field H c ^ 0.2 T. Upon 
H c , the (311), (400) magnetic Bragg peaks experience a 
continuous increase, showing a continuous polarization 
of Yb 3+ spin towards the direction of the magnetic field. 
The refinement of the 0.3 K pattern measured under 2 T 
actually yields a SF state with M=(±0.31(5), ±0.31(5), 
1.57(9 ))fiB in the global coordinate frame. The critical 
field is also confirmed by the ac magnetization measure¬ 
ment (Fig. 2(d)). At 75 mK, the ac susceptibility first 



FIG. 3: (color online) The electronic magnetic specific heat 
Cm for (a)Yb 2 Ge 207 and (b)Er 2 Ge 2 07 - The red dash lines 
show linear fits of the arrow-marked regions and the blue solid 
lines show fits considering the spin-wave gap. 


shows a peak at 0.12 T due to the domain alignment, 
and then another peak around H c = 0.22 T to enter 
the polarized state. With increasing temperature, both 
peaks’ positions move to lower fields and finally disappear 
above T Jy. This double peak feature is similar to that 
of EroGeoQTpS]. Along with our previous reported ac 
susceptibility data on Yb 2 Ge 2 07 | 29 j, a magnetic phase 
diagram is plotted in Fig. 2(e). However, due to the 
weak magnetic signal at (220) and the small H c , it’s dif¬ 
ficult to study how exactly this domain alignment affects 
the ( 220 ) peak’s intensity, which obstructs us to distin¬ 
guish between ^ 2 and ^ 3 . One noteworthy feature is that 
the dc magnetization measured at 0.6 K for Yb 2 Ge 2 07 
reaches 1.6 hb at 5 T. This value is consistent with that 
of Yb 2 Ti 2 07 and confirms the similar CEF scheme be¬ 
tween the Ge and Ti samples [32]. 

The selection of either ^2 or ^3 phase breaks the 
continuous U(l) symmetry, which requires a pseudo- 
Goldstone mode with a small spin-wave gap below T/y. 
For Er 2 Ti 2 07 , the inelastic neutron scattering has con¬ 
firmed the existence of this gap (~ 50 /ieV) |T8] . Mean¬ 
while, the specific heat data can reveal the information of 
this gap. Fig. 3(a) shows the electronic magnetic specific 
heat {Cm) ofYb 2 Ge 2 0 7 (details about obtaining C m are 
listed in the supplemental materials). Below T/y, C m fol¬ 
lows an almost prefect T 3 behavior down to 0.2 K, as 
the red dash line shows. However, it’s obvious that C m 
deviates from this straight T 3 line to a lower value be¬ 
low 0.2 K. Contrasting to a Goldstone mode where the 
Cm strictly follows a T 3 law, the gap that exists in the 
pseudo-Goldstone mode will multiply a component Ia{T) 
to T 3 , which is temperature dependent only in the tem¬ 
perature region that is comparable to the energy gap A. 
The relationship between the C m and A has already been 
derived in the supporting material of Ref. m- Here we 
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rewrite it as: 


ct 


MAkWa* ( 15 [°° X 2 (X 2 + (5 2 ) \ 3 

120^ 3 ^16tt 4 7 0 sinh 2 VS ) 

AI a (T)T 3 (1) 


where A/a is the Avogadro constant, ^ is the Boltz¬ 
mann constant, a is the lattice constant, v is the geo¬ 
metric mean of magnon velocity, X = /3k and 5 = (3A 
(dimensionless). The integration Ia{T) can be evaluated 
numerically with a given A. Ia(T) approaches a unity 
at high temperatures but decreases quickly when ksT is 
comparable to A, which leads the deviation of the C ^ 
from the T 3 behavior at low temperatures. The best fit 
of the measured to Eq. [I] with the A and A as two 
variables (blue line in Fig. 3(a)) yields the A = 24 jueV 
and A = 15.67 J.K -4 mol -1 , which corresponds to v = 
45.8 m/s. 

Similar analysis of the C m for Er 2 Ge 2 07 (Fig. 3(b)) 
yields a spin-wave gap A = 45 /ieV with A = 1.85 
J.K _4 mol _1 (corresponds to v = 132 m/s). One noticed 
feature is that at high temperatures, C m follows a T 2,72 
(not strict T 3 ) behavior. This could be due to the error 
bar introduced by the low temperature nuclear Schottky 
anomaly subtraction. 

With the decreasing lattice parameter or the increasing 
chemical pressure through the Sn to Ti to Ge samples, 
the magnetic ground states change accordingly (Table 
I). Given the fact that in these XY pyrochlores, the J ex 
dominate the magnetic properties, the chemical pressure 
can finely tune the J ex to lead to various magnetic ground 
states. This change of J ex is supported by the systematic 
changes of the Curie temperature and ordering temper¬ 
ature for XY-pyrochlores listed in Table I. Most strik¬ 
ingly, this is the first time to experimentally confirm an 
AFM ip 2 or 3 phase in Yb-pyrochlores despite the appar¬ 
ently different dominant exchange interactions between 
Yb and Er-pyrochlores(Ising-like J zz for Yb-pyrochlores 
and the XY-planar J± for Er-pyrochlores). This finding 
indicates there is general rules to unify the various mag¬ 
netic ground states of all effective spin-1/2 pyrochlores. 

Recent theoretical studies have made significant ef¬ 
forts to unify the magnetic ground states of the XY- 
pyrochlores. Wong et al. [24j have scaled the J ex by J± 
as three variables( J zz / J±, J z ±/ J±, J±±/J±) and calcu¬ 
lated a two dimensional magnetic phase diagram with the 
fixed ratio of J zz /J± , which contains continuous phase 
boundaries among the PC, SF, 7/2 and ips phases (the 
selection between 7/2 and ip 3 phases is through ObD). By 
comparing to the exchange interaction values obtained 
from the inelastic neutron scattering measurements, they 
successfully located the two Ti samples ( J zz /J± ~ -0.5, 
Jz±/J± ~ 0, J±±/J± ~ 1.0 for Er 2 Ti 2 07 and J zz /J± ~ 
3.0, J z ±/J± ~ -2.7, J±±/J± ~ 1.0 for Yb 2 Ti 2 07 ) in the 
7/2 and SF phase, respectively. Although we are short 
of knowledge of the exchange interaction values of other 



FIG. 4: (color online) Magnetic ground state phase diagrams 
for (a)Er 2 B 2 07 series and (b)Yb 2 B 2 07 series adopted from 
Ref. [24]. The dash areas are for just for the illustration pur¬ 
pose. The trends for the chemical pressure effects are shown 
as the direction of the arrows. 


XY-pyrochlores, here we located them in the J zz /J± = 
-0.5 and the J zz /J± = 3.0 phase diagrams adopted from 
Ref. [24]. This is based on two facts: (i) the phase dia¬ 
gram areas and boundaries are similar to each other over 
a wide range value of J zz / J±; (ii) the chemical pressure 
can finely tune but may not dramatically affect the ratio 
of Jzz/J±- As shown in Fig. 4, with increasing chemical 
pressure, two general trends are obvious: (i) the ground 
state moves downwards from PC state in Er 2 Sn 207 to 
7/2 in Er 2 Ti 2 07 and then ips in Er 2 Ge 2 07 for the Er- 
pyrochlores in the J zz /J± = -0.5 phase diagram; (ii) 
the grounds states move rightwards from the SF state 
of Yb 2 Ti 2 07 to the 7/2 or 7/3 region of Yb 2 Ge 207 in the 
Jzz/J± “ 3.0 phase diagram. 

These two trends can be successfully unified by the sce¬ 
nario that the increasing chemical pressure enhances J ±. 
For Er-pyrochlores with dominant XY type interactions, 
J zz and J z ± will take small values. Therefore, the in¬ 
creasing J± will primarily decrease the ratio of J±±/ J± 
to result in a downwards movement of the ground state. 
On the other hand, for Yb-pyrochlores with dominant 
local [111] Ising like interactions, J± and J±± will take 
small values. Therefore, the increasing J± will mainly 
decrease the ratio of J z ±/J± to result in a rightwards 
shift of the ground state to reach the AFM state for 
Yb 2 Ge 207 . Although without the values of the exchange 
interactions for all XY-pyrochlores, we cannot conclude 
the increase of the J± as the only reason for the change of 
ground states, the comparison between the reported J± 
values of Er 2 Sn 207 (J± = 1.35 meV) 22 and Er 2 Ti 207 
(J± = 6.7 meV) [14] supports our proposed scenario. 

Similar to Er 2 Ti 2 07 , the debate arises over what is the 
microscopic mechanism that breaks the continuous U(l) 
symmetry and selected the ordered phase below T/v- For 
the ObD scenario, the selection comes from the quan¬ 
tum fluctuations and is delicately tuned by the exchange 
parameters Jex [T4llT7] fFig. 4). For the CEF-induced 
energetic selection scenario considering additional CEF 
Hamiltonian and dipolar interaction ( Jdip ), the 7/2 phase 
is predicted in the Ts manifold pa Eg. The ip 3 phase 
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can only be achieved while adding a relatively strong 
Dzyaloshinskii-Moriya interaction ( Jdm ) EH]- Since the 
experimental results show that the change of lattice pa¬ 
rameter under chemical pressure has larger influences on 
J ex than CEF, Jd% v and Jdm [29]. The selection of dif¬ 
ferent ground state in E^TU 07 ( 7 / 2 ) and Er 2 Ge 207 (^ 3 ) 
seems to favor the ObD scenario. Furthermore, it is no¬ 
ticed that the values of magnon mean velocity and the 
gap in Yb 2 Ge 207 (v = 45.8 m/s, A = 24 /ieV) are both 
smaller than that of Er 2 Ge 207 (v = 132 m/s, A = 45 
/ieV), which is consistent with the ObD mechanism. A 
smaller v suggests a softer low lying mode in the spin 
wave spectrum that will result in a smaller energy dif¬ 
ference of spin-wave spectrum between the 7/2 and 7/3 
phases[14], for which a smaller gap value is expected. 
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Supplemental material 

1. EXPERIMENTAL SETUPS 

Polycrystalline sample Er 2 Ge 2 07 and Yb 2 Ge 207 were 
synthesized by the high-pressure and high-temperature 
(HPHT) technic. Proper ratio of starting materials 
Yb 2 03 , Er 2 03 , Ge 02 were mixed and synthesized in a 
Walker-type mutlianvil module (Rockland Research Co.) 
under 7 GPa and 1300 K. The ac susceptibility was mea¬ 
surements down to 20 mK on a home-made set up at 
SCM 1 of National High Magnetic Field Laboratory. The 
dc magnetization measurements were performed using 
a Quantum Design superconducting interference device 
(SQUID) magnetometer using a magnetic field of 0.01 T. 
The low temperature specific heat measurements were 
made in a Dilution Refrigerator option of the Physical 
Property Measurement System (PPMS, Quantum De¬ 
sign) using a standard semi-adiabatic heat pulse tech¬ 
nique. Elastic neutron scattering measurement was per¬ 
formed at Neutron Powder Diffractometer (HB-2A) at 
High Flux Isotope Reactor(HFIR) in Oak Ridge National 
Laboratory (ORNL). Neutron wavelength A = 2.41 A was 
used to maximize low angle magnetic scattering. Lat¬ 
tice and magnetic structure were refined though software 
package Fullprof-suite. 

2.MAGNETIC DOMAIN ALIGNMENT IN THE 
MAGNETIC FIELD 

For our powder neutron diffraction experiment on 
polycrystalline samples of Er 2 Ge 2 07 and Yb 2 Ge 2 07 , 
there are three main effects under consideration. (i)The 
grains are randomly oriented in the sample which deter¬ 
mine different orientations of the tetrahedrons. (ii)In the 
r 5 manifold, the U(l) rotational symmetry allows the 
four spins in a tetrahedron to rotate simultaneously in 
the local XY plane (we notate a as the angle of the spin 
relatively to the local x-axes in Fig. 1(d)). Then a single 
a defines an unique set of spin configuration. For a spe¬ 
cific 2 or V ^3 phase, there are six magnetic domains with 
a = n 7 r /3 (n = 0,...,5) for the ^2 phase and a = n 7 r /3 + 
7t/6 (n = 0,...,5) for the ^3 phase. In both cases, the six 
domains will be equally populated in a zero-field cooled 
sample below T/v- (iii)For a powder Bragg peak, it’s in¬ 
tensity is composed of all equivalent reflections. Specif¬ 
ically for the (220) Bragg peak, its intensity is equally 
contributed by six equivalent reflections: qi = 220 , q 2 
= 202 , q 3 = 022 , q 4 = 220 , q 5 = 202 , q 6 = 022 . 

We will begin with reflection of qi = 220 and show 
that the magnetic domain alignment behavior will be 
the same for all six reflections given a ^2 or ^3 phase. 
Namely, all reflections’ intensities will increase for a ^2 
phase and decrease for a ^3 phase. Thus the observed 


(220) Bragg peaks’ intensity drop in Er 2 Ge 2 07 under a 
small magnetic field identifies the ^3 phase. 

Fig. 1(a) illustrates our setup for the neutron powder 
diffraction experiments. The samples were pressed into 
pellets, wrapped by aluminum foils and fixed inside an 
aluminium can during the experiment to prevent mechan¬ 
ical motion of the grains. The magnetic field was applied 
vertically so that it is perpendicular to the scattering 
plane. For the qi = 220 reflection, its intensity comes 
from all grains with its [110] direction along qi. Then 
the geometry defines a grain-average plane that is per¬ 
pendicular to qi where the magnetic filed direction, all 
grains’ [001] and [110] axes he within (Fig. 1(b)). Assum¬ 
ing equally distribution of different grains in the powder 
sample, then the effect of averaging different grains with a 
unique magnetic field direction in the grain-average plane 
(the case of the experiment) is equivalent to that of aver¬ 
aging random magnetic field directions to a single grain 
in a field-average plane (the case for modeling). The re¬ 
sulting (hhl) field-average plane following the definition 
of the grain orientation is shown in Fig. 1(c). Its re¬ 
lationship to the Er 3+ tetrahedron is illustrated in Fig. 
i(d). 

Fig. 2 shows the magnetic Bragg peaks’ intensity as 
a function of a for all (220) type reflections calculated 
using the software package Fullprof-suite. For the qi = 
220 reflection, the intensity for the six -02 domains (dots 
in Fig. 3(a)) will take two large values (for a = 0 and i r 
)and four small values (for a = 7r/3, 27 t/3, 4tt/3, 57t/ 3). 
With the magnetic field plane defined above for qi, the 
a = 0 and a = tt states will be selected by the symmetry 
as shown in the first row of Fig. 3. Then the intensity 
of qi = 220 reflection will experience an increase due to 
the magnetic domain alignment under magnetic fields. 
On the other hand, the intensity for the six -03 domains 
(square in Fig. 2(a)) will take four large values (for a = 
7r/6, 57 t/6, 77t/6, 117t/ 6) and two small values (for a = 
7t/2, 37t/ 2). The later two spin configurations with small 
intensities will be selected by symmetry as shown in the 
bottom row of Fig. 3, resulting the intensity drop for the 
qi = 220 reflection. 

The same analysis will apply for other reflections qi 
(i = 2,...,5) which define different magnetic planes (Fig. 
l(e-i)) with according magnetic reflection intensities (Fig. 
2(b-f)). Take the q6 = 022 for example, the according 
(lhh) field-average plane (Fig. l(i)) favors <a = 7r/3, 47 t/ 3 
states for the ^2 phase and a = 57r/6, 117r/6 states for 
the -03 phase. The resulting intensity change in Fig. 2(f) 
will be exactly same as the qi situation. 

In summary, considering the grain population and 
equivalent reflections effect in the powder sample, the 
decrease of (220) Bragg peak’s intensity due to magnetic 
domain alignment identifies the ^3 phase in Er 2 Ge 2 07 . 
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FIG. 5: (a)Experimental setup for the powder neutron diffraction experiment. The six equivalent (220) type reflections are 

labeled as qi (i = 1,...,5). (b)(c)Illustration of the grain-average and the according field-average for a single reflection qi. 
(d-i)Different magnetic field-average planes for different qi. 
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FIG. 6: Magnetic Bragg peak’s intensity as a function of ro¬ 
tation angle a for different (220) type reflections. The solids 
dots represent the six ^2 domains and the open squares rep¬ 
resent the six ip 3 domains. 
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FIG. 7: The six magnetic domains’ spin configurations for 
the and ^3 phases. The magnetic field-average plane is 
illustrated for qi = 220. 


3.ELECTRONIC MAGNETIC CONTRIBUTION 
TO THE SPECIFIC HEAT 


We obtain the electronic magnetic contribution (C m ) 
to the specific heat for Er 2 Ge 207 and Yb 2 Ge 2 07 after 
subtracting the lattice (Cl) and magnetic nuclear (Cat) 
contributions from the total measured specific heat (C t) 
of each sample. The lattice contribution was estimated 
measuring the structurally similar non-magnetic material 
Lu 2 Ge 2 07 (see Fig. 4). The exact expression for the 


nuclear specific heat is [ll: 


Cn = 
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E E {wf-WiWjy 


( k B Ty 


i — I j=I Wi + Wj 

E E 

i=-ij=-i 


( 2 ) 



























where I is the nuclear spin and the energy levels Wi are 
given by pQ: 
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Here a’ = /&#/ and P = 3e 2 Qq/4k/(2/-l) are the 

magnetic interaction parameter and the quadrupole cou¬ 
pling constant. We use the nuclear spin I of each com¬ 
pound and fit the low temperature part of the data with 
a’ and P as adjustable parameters to obtain the nuclear 
contribution to the specific heat as shown in Fig. 4(a)(b). 

The C m/ T and the integrated magnetic entropy for 
both samples are plotted in Fig. 4(c)(d). For both sam¬ 
ples, the recovered magnetic entropies almost reach Rln2. 


FIG. 8: Total specific heat (Ct) for (a)Yb2Ge207 and (b) 
Er2Ge207- The electronic magnetic contribution (C m ) is ob¬ 
tained after the subtraction of the lattice (Cl) and magnetic 
nuclear (Cjv) specific heat. C m /T and the integrated entropy 
for (c)Yb2Ge2 07 and (d)Er2Ge207- 


[1] O. V. Lounasmaa, A. J. Freeman e R. B. Frankel, ”Nuclear 
specific heats in metals and alloys,” in Hyperfine Interac¬ 
tions , New York, Academic Press, 1967, pp. 467-496. 





























